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Abstract. The abundances of O, Mg, Si and Fe in the intracluster medium of the Centaurus cluster are derived. 
The Fe abundance has a negative radial gradient. In solar units, the Si abundance is close to the Fe abundance, 
while the O and Mg abundances are much smaller. The high Fe/O and Si/O ratio indicate that metal supply from 
supernovae la is important and supernovae la synthesize Si as well as Fe. Within 2', the O and Mg abundances are 
consistent with the stellar metallicity of the cD galaxy derived from the Mg2 index. This result indicates that the 
central gas is dominated by the gas from the cD galaxy. The observed abundance pattern of the Centaurus cluster 
resembles to those observed in center of other clusters and groups of galaxies. However, the central Fe abundance 
and the Si/Fe ratio are 40 % higher and 30% smaller than those of M 87, respectively. Since the accumulation 
timescale of the supernovae la is higher in the Centaurus cluster, these differences imply a time dependence of 
nucleosynthesis by supernovae la. 
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1. Introduction 

The intracluster medium (ICM) contains a large amount 

■ of metals, which are synthesized by supernovae (SNe) la 
and SNe II in cluster galaxies. A knowledge of the con- 
tributions from the two types of SNe to the metal abun- 
dances are important to understand the origin of the met- 

■ als and to study the evolution of galaxies in clusters. 

The XMM-Newton observatory enables us to derive 
O and Mg abundances, which are not synthesized by SN 
la. For the ICM of M 87 observed with the XMM, the 
O/Si ratio is found to be about a half of the solar ratio, 
while Fe and Si have similar values in units of the solar 
abundance (Bohringer et al. 2001; Finoguenov et al. 2002; 
Matsushita et al. 2003). Such low O abundances compared 
to the Si abundances are also derived for the centers of 
other clusters and groups of galaxies (e.g. Tamura et al. 
2001b ; Boute et al. 2003, and others). 

In order to account for the observed O/Si/Fe abun- 
dance pattern of M 87, the Fe/Si ratio of ejecta of SN la 
should be about 1 solar ratio. This value is significantly 
higher than expected and probably related with dimmer 
SNe la observed in old stellar systems (Hamuy et al. 1996; 
Ivanov et al. 2000). Since the mass of synthesized Ni 56 de- 
Send offprint requests to: K. Matsushita 
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termines the luminosity of each SN, the abundance ratio 
in the ejecta may be correlated to the luminosity of SN 
la. In this context a lower Ni 56 yield implies less complete 
nuclear burning with a higher ratio of a-elements such as 
Si to Fe. The observed O/Si/Fe pattern of the Galactic 
stars also indicates that the Si/Fe ratio in ejecta of SN 
la depends on stellar metallicity, i.e., age of the system 
(Matsushita et al. 2003) 

The variation in the observed abundance pattern re- 
flects variations in the explosion models, e.g. a high 
fraction of Fe for the classical deflagration model, W7 
(Nomoto et al. 1984) or a larger fraction of Si for the 
delayed detonation model, WDDs (Nomoto et al. 1997; 
Iwamoto et al. 1999). The W7 model, which can well re- 
produce observed optical spectra of SNe la (Nugent et al. 
1997), predicts a high Fe/Si ratio of 2.6 in units of the so- 
lar ratio. In contrast, the WDD models give a wider range 
of Fe/Si ratios from 1 to 3 in units of the solar ratio, which 
might be related to the age of the system (Umeda et al. 
1999). 

In this paper, we report on the O/Mg/Si/Fe abun- 
dance pattern of the Centaurus cluster observed with 
XMM. The temperature properties of the cluster are de- 
scribed in Takahashi (2004). The Centaurus cluster is a 
nearby cluster, centered on the cD galaxy, NGC 4696. 
A positive temperature gradient and a negative abun- 
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dance gradient are observed with ASCA (Fukazawa et 
al.1994; Ikebe et al.1999). The cluster is expected to have 
a cooling flow (Allen & Fabian 1994), and Chandra ob- 
served an unusual plume-like cool structure at the center 
(Sanders & Fabian 2002). However, as in other cooling 
flow clusters, the luminosity of the central cool compo- 
nent is much smaller than the expected value from the 
cooling flow model (e.g. Makishima et al. 2001; Tamura et 
al. 2001a; Bohringer et al. 2002; Matsushita et al. 2002). 
The problem of the missing of the cooling component be- 
comes an advantage in determination of elemental abun- 
dances, since the observed narrow temperature range of 
the ICM strongly reduces the ambiguity in the abundance 
measurement. 

We adopt for the solar abundances the values given by 
Feldman (1992), where the solar Fe abundance relative to 
H is 3.24 xl0~ 5 in number. A value for the Hubble con- 
stant of Hq = 70 km/s/Mpc is assumed. Unless otherwise 
specified, we use 90% confidence error regions. 

2. Observation 

The Centaurus cluster was observed with XMM-Newton 
on January 3rd, 2002. The thin filter is used for the 
EMOS and the EPN. We have selected events with pat- 
tern smaller than 5 and 13 for the EPN and the EMOS, 
respectively. The background spectrum was calculated for 
each spectrum by integrating blank sky data in the same 
detector regions. Among deep sky observations collected 
with XMM, we have selected the data with the most sim- 
ilar background to that of the Centaurus cluster, after 
screening background flare events in the data and the 
background following Katayama et al. (2004). The effec- 
tive exposures of the EPN and the EMOS are 29 ks and 
32 ks, respectively. The spectra are accumulated within 
rings, centered on the X-ray peak of the cluster. After 
subtracting the background and correcting for the effect 
of vignetting, deprojected spectra are calculated assuming 
spherical symmetry. For the EPN spectra, the response 
matrix file and the auxiliary response file corresponding 
to each spectrum are calculated using SAS-v5.4.1. For the 
EMOS data, we used ml_thinlv9q20t5r6_all_15.rsp. The 
spectral analysis uses the XSPEC_vl 1.2.0 package. For the 
EPN data, we used the spectral range from 0.5 to 7.3 keV, 
in order to discard strong instrumental lines around 8 keV. 
We fitted the EMOS data in the spectral range from 0.4 
to 9.0 keV. At r > 8', the energy range of 1.45-1.55 keV 
of the EMOS is discarded due to a strong instrumental 
line. Here, r is the projected radius from the X-ray peak. 
We also denote R as the deprojected or three-dimensional 
radius. 

3. Spectral fit with the MEKAL model 

The first step of the analysis is fitting the annular (pro- 
jected) spectra with a single temperature MEKAL (e.g. 
Kaastra 1992; Liedahl et al. 1995) model, and the depro- 
jected spectra with a single and two-temperature MEKAL 
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Fig. 1. The deprojected spectrum at i?=0.5-2' of the 
EMOS1+EMOS2 of the Centaurus cluster fitted with the 
single temperature MEKAL model (black thin lines), the 
two-temperature MEKAL model (black dotted lines), the 
multi-temperature MEKAL model (black bold lines), and 
the multi-temperature APEC model (gray bold lines). 



model with photoelectric absorption. Abundances of C 
and N are fixed to be 1 solar and those of other elements 
are determined separately. The abundances of each ele- 
ment for all temperature components are assumed to have 
the same values. In addition, the deprojected spectra are 
fitted with a multi-temperature model which is a sum of 
14 temperature components. Since a spectrum from an 
isothermal plasma can be well reproduced by a sum of 
spectra of plasma of two neighboring temperatures with 
similar elemental abundances, we selected temperatures in 
order to derive metal abundances with an accuracy within 
5 %. Thus, in order to determine the abundance distribu- 
tion, this multi-temperature model is enough to fit any 
temperature distribution at least when all temperature 
components have the same abundances. 

Table 1 and 2 summarize the results. Even for the 
deprojected spectra, the single temperature model is 
not acceptable at R < 2'. The two-temperature model 
dramatically reduces the x 2 , and the multi-temperature 
model slightly improves x 2 - Outside this radius, x 2 de- 
rived from the two- and the multi-temperature model 
are slightly smaller than those of the single-temperature 
model. Figure 1 shows a representative spectrum for 
i?=0.5-2'. Using the two- and the multi-temperature 
MEKAL model, not only the continuum and the Fe-L 
bump, but also the strength of the Fe-K line is fitted 
well. The consistency between the Fe-L and the Fe-K will 
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Table 1. Results of spectrum fitting of the projected spectra. 
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0.05- 0.25 
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0.00 
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EPN 
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EPN 
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be discussed in detail at H4.3.3I There remains a residual 
structure around 1.2—1.5 keV at the Fe-L/Mg-K region. 
Therefore, in Table 1 and 2, we do not show errors in the 
Mg abundances. 

The temperature of the ICM decreases sharply toward 
the center (Figure |2J , and we need more than two tem- 
perature components to fit the deprojected spectra at the 
central region. The observed sharp temperature gradient 
means that the single temperature model is not enough to 
fit each deprojected spectrum, even if the ICM is isother- 
mal locally. In addition, the cluster is not circularly sym- 
metric and the cool plume-like structure at the center 
(Sanders & Fabian 2002) might introduce an additional 
temperature component. 

The hydrogen column density, Nh is consistent with 
the Galactic value at R > 1', although the EPN and 
the EMOS give slightly different results by a few times 
10 20 cm~ 2 (Figure Within R < 1', the two- and the 
multi-temperature models give similar values to those of 
R > V . Higher column densities derived from the single 
temperature model fit within R < V are not acceptable 
due to their high \ 2 ■ Therefore, Nh is consistent with the 
Galactic value within the whole field of view. 

Discarding the results of unacceptable fits (reduced 
X 2 > 2), the Fe abundance is ~0.6 solar at R > 8' 
and within the radius, it increases steeply toward cen- 
ter. Within 2', the gradient becomes weaker. The average 
value of the Fe abundance at R < 2' is 2.3 solar (Figure|2I). 
The EPN and the EMOS give consistent Fe abundances 



within 10% (Table 2). The Fe abundances from the de- 
projected spectra are slightly higher than those from the 
projected spectra. 

At each radius, the Si abundance is close to the Fe 
abundance (Figure [21 . The uncertainties in the Si/Fe 
ratio are better constrained than those in the Si and Fe 
abundances themselves (@J, since the derived abundance 
of Si and Fe are correlated to each other. On average, the 
Si/Fe ratio is 0.89±0.04 and 0.75±0.04 in units of the solar 
ratio using the EMOS and EPN, respectively. 

The O abundances are much smaller than the Si and Fe 
abundances (Figure|21 4) . The two detectors give inconsis- 
tent O abundances. Using the EMOS, the O abundance is 
about 0.8 solar at R—0. 5-2.0' and at R >2', it reduces to 
0.4 solar (Figure EJ), and the O/Fe ratio increases toward 
outer regions (Figure In contrast, the EPN gives a sig- 
nificantly smaller O abundance than the EMOS especially 
in the outer regions (Figure|3Jl, and the profile of the O/Fe 
ratio is consistent with no radial gradient (Figure |3J. 

In summary, as in center of other clusters and groups of 
galaxies, the Fe and Si abundances are close to each other 
with strong negative radial gradients. The O abundance 
is at least a factor of 2 smaller. The discrepancies in the 
O and Si abundances between the EMOS and EPN will 
be discussed in section 4.3. 
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Table 2. Results of spectrum fitting of the deprojected spectra 
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IT 
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7.7 
7.8 
8.3 



0.5 
+0.6 
0.5 
+0.9 
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4. Uncertainties in the abundance determination 

In the abundance determination, we must examine uncer- 
tainties such as the dependence on the plasma code and 
on the background subtraction. In this section, we study 
these uncertainties and also calibrate the abundances with 
the strengths of the lines. 

4.1. The effect of the plasma code 

In order to constrain the effect of the plasma code on the 
abundance determination, we compared the results from 
the APEC (Smith et al. 2001) model with those from the 
MEKAL model. We used the APEC version 1.1 model 
in XSPEC version 11.2. As in the MEKAL model fit, we 
fitted the deprojected spectra with a single, and a multi- 
temperature APEC model. 

The results are summarized in Table 3. With the 
APEC model, the fit to the Fe-L/Mg-K region at 1.2-1.5 



keV improved from that with the MEKAL code (Figure 
1). As a result, the reduced \ 2 are consistent or slightly 
smaller than those from the MEKAL model fit. The de- 
rived temperatures, hydrogen column densities, and the Si 
and Fe abundances are consistent within 10% with those 
from the MEKAL code (Figure |3 EJ) . 

However, the O and the Mg abundances increased by 
-20% and- 80%, respectively, from the MEKAL model 
fit. The difference in the O abundance is consistent with 
the difference of the strength of the H-like O line between 
the APEC and MEKAL models, since at 2 keV, that of 
APEC is 10% larger than that of MEKAL model, and at 4 
keV, the difference increases to 30%. The difference in the 
Mg abundance is not due to the difference in the strength 
of the Mg-K line between the two codes which is smaller 
than that of O below 4 keV. It should be caused by the 
difference in the Fe-L modeling around the Mg-K line. We 
will discuss it later in §4.3.2. 
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Fig. 2. The radial profiles of the temperature and the hydrogen column density. Dotted and solid lines correspond 
to the results of the projected spectra and the deprojected spectra, respectively. Black, red and green correspond to 
the single, two, and the multi-temperature MEKAL model, respectively. The results from the EMOS are shown in 
diamonds and those from the EPN are shown in crosses. The dashed horizontal line in the profile of the hydrogen 
column density represents the Galactic value. 

Table 3. Results of spectrum fitting of the deprojected spectra within 8' and projected spectrum at 8-13.5' of the 
EMOS with the single- and the multi-temperature APEC model. 
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4.2. Effect of the background uncertainty 

One of the most severe problems in analyzing XMM- 
Newton data are uncertainties in the background. In this 
paper, we adopted deep survey data as a background, 
which agree well with the data of the Centaurus cluster 
above 8 keV (Figure 0). In addition, in order to study 
the effect of the uncertainty of the background, we artifi- 
cially scaled the background by 10%, and fitted the pro- 
jected EMOS spectra at r > 2' in the same way as previ- 
ously. Since spectral variations of the background amount 
to about 10% below 1 keV when the background is scaled 
with the count-rate outside the field of the view of the 
detector (Katayama et al. 2004), and the relatively high 
Galactic A^h towards the Centaurus cluster absorbs the 
low energy photons from the Cosmic X-ray Background 
(CXB). 

The results are summarized in Figure |S1 Within 8', 
the derived temperatures, the Fe abundances, the O/Fe 
and the Si/Fe ratios are consistent within a few % with 



those derived from the background used in §3.1. Outside 
8', the 10% change of the background changes derived 
parameters by ~ 10%. Therefore, the background problem 
is not severe in this observation. 

4.3. Comparisons between the line strength and model 
fit 

4.3.1. The line strengths of O and Si and the 
abundance determination 

In order to study the discrepancy in the O and the Si abun- 
dances between the EMOS and the EPN data, it is impor- 
tant to check their line strengths directly. Complications 
arise since the energy resolution of the EMOS has be- 
come worse due to the radiation damage in orbit (Kirsh 
2002), and that of the EPN depends on the distance from 
the read-out. In addition, there are still discrepancies of 
several % in the response matrix between the EPN and 
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Fig. 3. Radial profiles of the abundances of Fe, Si and O. 
The meanings of symbols are the same as those of Figure 
2. The data with a reduced x 2 > 2 are not shown. 



EMOS, which might give a large uncertainty in strengths 
of weak lines. 

The strengths of Ka lines of H-like O and Si are de- 
rived from fittings of the spectra around the lines with 
Gaussians. The details are the same as in Matsushita et 
al. (2003). First, the line widths of the Gaussians were 




R (arcmin) 

Fig. 4. Radial profiles of the O/Fe (black lines) and the 
Si/Fe (gray lines) ratios in units of the solar elemental 
ratio derived from the multi-temperature MEKAL model 
of the EMOS (diamonds) and the EPN (crosses) spectra. 
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Fig. 5. The Fe abundance profile derived from the 
MEKAL (gray data) and the APEC (black data) model, 
from the single-temperature model (crosses) and the 
multi-temperature model (diamonds). The data with a re- 
duced x 2 > 2 are not shown. 



frozen to zero, and then allowed to be free. Figure [5] com- 
pares the best fit MEKAL model and the Gaussian fits 
for representative spectra at i?=2-8'. In the EMOS spec- 
trum, we can see a weak bump of the H-like O line, and a 
clear bump of the H-like Si line. In contrast, in the spec- 
trum of the EPN, which has an energy resolution worse 
than that of the EMOS, the continuum and the O line are 
hardly distinguished and the Si line is much broader than 
the EMOS one. 

For the EMOS spectra, the line strengths of K Q of H- 
like O and Si from the Gaussian fits are consistent within 
10% with the those from the MEKAL model, except for 
the O line strength within 0.5', where the O line is com- 
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Fig. 6. The ratio of the 0(open triangles), Mg(crosses), 
Si(open circles) and Fe(closed circles) abundances and the 
reduced x 2 (closed boxes) between the multi-temperature 
APEC and the multi-temperature MEKAL model. 
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Fig. 7. The raw EMOS spectrum at r > 8' (black crosses) 
and the adopted background spectrum (black diamonds), 
and the background subtracted spectrum (gray crosses). 

pletely hidden in the Fe-L lines (Figure ll(J|) . This result 
indicates that the effect of the decrease of the energy reso- 
lution of the EMOS detector should not affect the derived 
O and Si abundances. 

In contrast, allowing the line width to be a free pa- 
rameter, the Gaussian fits to the EPN spectra give sys- 
tematically larger line strengths than the MEKAL model 
fit. As a result, adopting the values from the Gaussian fits, 
the EPN abundances become consistent with the EMOS 
ones. Therefore, the O and Si abundances derived from the 
EMOS fit should be reliable, since it has a better spectral 
resolution and it gives consistent strengths of the Ka lines 
between the Gaussian fits and the MEKAL model fits. 

As shown in the Figure ^2 the ratio of strengths of 
Ka lines of H-like O and Si becomes constant above 1.7 
keV and increases toward lower energies. This temper- 
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Fig. 8. The ratio of derived parameters of the projected 
EMOS spectra with the multi-temperature model using a 
background scaled by a factor of 1.1 to those using the 
unsealed background used in the paper. 
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Fig. 9. Deprojected spectra of the EMOS 

(EMOS1+EMOS2; lower curve) and the EPN (up- 
per curve) at R=2. 0-8.0' around the energy of the lines of 
O (left panel) and Si (right panel). The gray data points 
and black solid lines correspond to the best fit model with 
one temperature MEKAL model and that with Gaussians 
with line widths, respectively. The black dotted lines 
represent the continuum emissions as described by the 
Gaussian fit. 

ature dependence means that any temperature distribu- 
tion cannot yield a O/Si ratio larger than 0.6 solar within 
4'. Furthermore, we got a higher value of the line ratio 
at R > 8' by a factor of 2. It implies that the O/Si ra- 
tio increases in the region, since the contribution to the 
O line from the component below 1.7 keV of the multi- 
temperature MEKAL model is smaller than 10%. 

4.3.2. The line strength of the Mg-K 

It is also important to derive the strengths of the Ka-lines 
of Mg directly, since the Mg abundances derived from the 
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Fig. 10. The ratio of line strengths Ka of H-like (upper 
panel) and Si (lower panel) derived from the Gaussian fit 
to those of best fit values from spectral fit for the EMOS 
(black lines) and EPN (gray lines). Crosses and diamonds 
correspond to the fitting with a Gaussian with line width 
and free line width, respectively. Errors correspond the 
68% level. 
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Fig. 11. The line ratios of Ka lines of H-like O and Si 
of the projected spectra of the EMOS, plotted against 
the best fit temperature fitted with the single MEKAL 
model. Dashed lines correspond to constant abundance 
ratios in units of the solar ratio. Errors represent the 68% 
confidence level. 

spectral fits depend on the modeling of the underlying 
Fe-L lines. We therefore fitted the EMOS spectra within 
the energy range of 1.3-1.65 keV with Gaussians plus one- 
two- or three-temperature MEKAL or APEC model. 

All of the Gaussian models can fit the spectra around 
the Mg-K/Fe-L region well (Figure H2J. The derived 
strengths of the Ka-line of H-like and He-like Mg do 
not depend on the temperature structure (i.e. one-, two, 
or three-temperature model) or the plasma code (i.e. 
MEKAL or APEC model). The line strengths agree well 
with those derived from the fitting the spectra described in 
§4.1, but are systematically higher than those derived from 
the multi-temperature MEKAL model described in §3, 




channel energy (keV) 

Fig. 12. The deprojected spectrum of the EMOS at 
R=0. 5-2.0' (black closed circles) fitted with the multi- 
temperature MEKAL model (green closed circles), the 
multi-temperature APEC model (red closed circles), the 
Gaussians with MEKAL models (green dashed lines) and 
the Gaussians with the APEC models (red dashed lines). 
The contributions of the continuum plus the Fe-L are also 
shown for the multi-temperature APEC model (red open 
circles) and the multi-temperature MEKAL model (green 
open circles) and the Gaussian fit (red and green dashed 
lines). 



which cannot reproduce the spectra around the Mg-K/Fe- 
L region (Figure EU- Therefore, it is reasonable to adopt 
the Mg abundances derived from the multi-temperature 
APEC model fit as the Mg abundances of the ICM. 

We note that the Mg/O ratio does not depend on 
the temperature structure since the line ratio of Ka lines 
of these two elements is nearly constant above 1.1 keV 
(Matsushita et al. 2003). 

4.3.3. The line strength of Fe-K 

Since the temperature dependence of the intensities of the 
Fe-L and Fe-K lines is completely different, a comparison 
between the Fe-L and the Fe-K lines strongly constrains 
the temperature structure and the Fe abundance. In the 
spectral fitting by the MEKAL model or the APEC model, 
the strengths of the 6.7 keV Fe-K lines are fitted well, 
although the fit to the Fe-L lines mostly determines the 
Fe abundance due to its high statistics (Figure ITlT) . 

When the Fe abundance is shifted by 40% in the fits 
with the multi-temperature model, the Fe-L structure can 
be roughly fitted although the x 2 increased significantly 
(Figure . However, the strengths of the Fe-K line can- 
not be reproduced. Therefore, the uncertainties in the 
Fe-L line fit, due to the temperature structure and the 
plasma code should not be severe problems at least in the 
Centaurus cluster. 

Since abundances derived from the Fe-K line may have 
an uncertainty due to the uncertainty of the background, 
we have also derived the line strength of the Fe-K line di- 
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Fig. 13. The profile of the line strengths of the Ka-line 
of He-like and H-like Mg of the best fit model of the 
multi-temperature model (horizontal solid lines) and de- 
rived from the Gaussians plus a single temperature model 
(solid diamonds), two temperature model (dashed dia- 
monds), and three-temperature model (dotted diamonds). 
The gray and black data correspond to the MEKAL and 
APEC model, respectively. Errors correspond the 68% 
confidence level. 
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Fig. 14. The \ distributions of the deprojected spectrum 
at #=0.5-2.0' fitted with the multi-temperature APEC 
model. The solid black line corresponds to the best fit 
model, and the dotted black line and dashed gray lines 
correspond to the model where the Fe abundance is fixed 
to be 1.6 solar and 3.2 solar, respectively. 



rectly from a Gaussian fit. The width of the Gaussian was 
allowed to be a free parameter, since the 6.7 keV Fe-K con- 
sists of several strong lines. In figure 1X31 the derived line 
strengths from the Gaussian fit are compared to the line 
strengths of the Fe-K line derived from the best fit multi- 
temperature MEKAL model. These values agree within 
10%. Thus, the line strengths of the Fe-K line are consis- 
tent with the Fe abundances derived from the MEKAL 
model fit and a change of the Fe abundance cannot repro- 
duce the observed line strength of the Fe-K. 
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Fig. 15. The ratio of line strengths of the 6.7 keV Fe-K 
line derived from the Gaussian fit to those of best fit values 
from spectral fit for the EMOS (black lines) and EPN 
(gray lines), where the Fe-L is important to determine the 
Fe abundance. Errors correspond to the 68% level. 

5. Discussion 

5.1. Summary of the results 

The abundances of O, Mg, Si, and Fe in the ICM of the 
Centaurus cluster are derived. We have checked uncertain- 
ties in the temperature structure, in the background sub- 
traction, and in the spectral model and studied their ef- 
fects on the abundance determination. The Fe abundance 
is ~ 2.3 solar within 2' and decreases to ~0.6 solar at 
R > 8'. The Si abundance is close to the Fe abundance. 
In contrast, the O and Mg abundances are much smaller 
than the Fe and Si abundances in units of the solar ratio. 

The abundance pattern of the ICM of the Centaurus 
cluster, the low O abundance compared to those of Si and 
Fe, resembles to those observed in center of other clusters 
and groups of galaxies, such as M 87 (Matsushita et al. 
2003), A496 (Tamura et al. 2001b), Fornax cluster (Buote 
et al. 2002) and group-center elliptical galaxies NGC 4636 
(Xu et al. 2002), NGC 5044 (Tamura et al. 2003; Boute 
et al. 2003). However, there are small differences in the 
observed abundance pattern of the ICM. For example, the 
central Fe abundance of the Centaurus cluster (2.3 solar) 
is 40% larger than that of M 87 (1.6 solar). In addition, the 
Si/Fc ratio of the former is 20% smaller than that of the 
latter. In this section, we mainly compare and discuss the 
abundance pattern of the Centaurus clusters with that of 
M 87, since these two clusters are the nearest and the best 
studied objects among relatively low temperature clusters. 

5.2. Comparison with stellar metallicity and SN II 
abundance pattern 

The O and Mg in the gas around a cD galaxy is a mix- 
ture of those in the ICM and those in the accumulation of 
stellar mass loss from the cD galaxy, since these elements 
are only synthesized by SN II. The O and Mg abundances 
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Fig. 16. Mg (closed circles) and O (diamonds) abun- 
dances derived for the Centaurus cluster (black) and M 
87 (gray) from the multi-temperature MEKAL (dotted 
lines) and APEC (solid lines) model. The dashed lines rep- 
resent the stellar metallicity of NGC 4696 (black) and M 
87 (gray) derived from the optical Mg2 index (Kobayashi 
& Arimoto 1999). 



of the mass loss reflect those of mass losing stars, which 
contribute significantly to the optical light from elliptical 
galaxies. The stellar metallicity of elliptical galaxies are 
usually derived from the Mg2 index in the optical spec- 
tra, which depends mainly on the Mg abundance and stel- 
lar metallicity where O contributes most, and weakly on 
the age population of stars. Therefore, comparing the O 
and Mg abundances of the gas and stellar metallicity de- 
rived from the optical Mg 2 index is important to study the 
gas flow and contribution of stellar mass loss, and stellar 
metallicity of the cD galaxies. 

The Mg abundances of the two clusters are close 
to their stellar metallicity derived from the Mg2 index 
(Figure El Kobayashi & Arimoto 1999). The gradient of 
the Mg abundance of the gas around M 87 is steeper than 
in the Centaurus cluster. At a radius of 10-20 kpc, the Mg 
abundance of the Centaurus cluster is a factor of 2 larger 
than that of M 87 (Figure flT)|) . The stellar metallicity gra- 
dient of M 87 is also steeper than that of NGC 4696, the 
cD galaxy of the Centaurus cluster. The agreements of be- 
tween the Mg abundances with stellar metallicity at the 
central regions indicates that in these regions, the gas is 
dominated by gas ejected from the cD galaxies. 

The Mg/O ratio of the gas of the Centaurus cluster is 
1.2 ± 0.3 of the solar ratio on average as determined with 
the multi-temperature APEC model. This value is close to 
that of the ICM within 4' of M 87, which is 1.3± 0.1 solar 
ratio (Matsushita et al. 2003), and also to the value of 
1.3±0.2 solar ratio of NGC 4636 (Xu et al. 2003) derived 
from the RGS. The observed Mg/O ratios of the gas at the 
cluster centers are close to [Mg/O] of the Galactic stars 
(Edvardsson et al. 1993). Therefore, at least for the Mg/O 
ratio, the Galactic SN II and those trapped in stars in the 
cD galaxies have no difference. 
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Fig. 17. The Fe abundance of the ICM of the center of 
Centaurus cluster (black) and M 87 (gray). 



5.3. The Fe abundance profile and SN la contribution 

The Fe in the gas at the core of clusters is also a mixture 
of those in the ICM and the recent supply from the cD 
galaxy. The latter contains Fe synthesized by SN la and 
those coming from stars through stellar mass loss, since 
Fe is synthesized by both SN la and SN II. 

The Fe abundance of the ICM of the Centaurus cluster 
is systematically higher than that around M 87 at a same 
distance from the cluster center (Figure I17|l . The peak 
Fe abundance of the Centaurus cluster is 2.3 solar. This 
value is significantly higher than that of M 87, which is 
1.6 solar. Assuming that the O/Fe ratio of SN II is 3, the 
central Fe abundances from SN la of the Centaurus cluster 
and M 87 become 2 and 1.3 solar, respectively. Since the 
SN la contribution to the Fe abundance is dominant, the 
uncertainty in the assumption for the O/Fe ratio does not 
change the result very much. 

The Fe abundance of gas from SN la in an elliptical 
galaxy is proportional to M|j|0sn / 'a* (see Matsushita et 
al. 2003 for details). Here, M|n is the mass of Fe synthe- 
sized by one SN la, #sn is SN la rate, and a* is stellar 
mass loss rate. The higher Fe abundance of the gas in the 
central region of the Centaurus cluster indicates a higher 
MgN^SN/ a *- The cooling flow rate derived from the stan- 
dard cooling flow model is higher in the Centaurus cluster 
(Allen & Fabian 1994) than that of the M 87 (Matsushita 
et al. 2002). Bohringer et al. (2004) questioned the stan- 
dard cooling flow model comparing the Fe mass profile of 
4 clusters of galaxies with cooling cores including M 87 
and the Centaurus cluster. They concluded that we need 
long enrichment times (> 5 Gyr) in order to accumulate 
the observed Fe abundance peak even if the SN la rate 
was larger in the past. 

As shown in Figure 18, a difference between M 87 and 
the Centaurus cluster is the ratio of the Fe mass to the 
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Fig. 18. The Fe abundance (the upper panel) and 
(Fe/Si) SNIa (the middle and the bottom panel) in units 
of the solar values are plotted against the ratio of the 
Fe mass to B-band luminosity of the cD galaxy (IMLR), 
for M 87 (gray circles) and the Centaurus cluster (black 
circles), when adopting the abundance pattern of SN II 
from the average value of Galactic metal poor stars (the 
middle panel; Clementini et al. 1999) and from nucleosyn- 
thesis model (the bottom panel) by Nomoto et al. (1997) 
assuming Salpeter's initial mass function (Iwamoto et al. 
1999). 



5.4. Abundance pattern of SN la 

Si is both synthesized by SN la and SN II, and the abun- 
dance pattern of O, Si, and Fe strongly constrains the 
nucleosynthesis of the SN la when Fe is dominated by SN 
la. Figure ITH1 summarizes the ratio of these three elements 
in the ICM of the Centaurus cluster and M 87 (Matsushita 
et al. 2003). As in M 87, the data of the Centaurus cluster 
cannot be explained by a sum of SN II abundance pattern 
in the Galactic metal poor stars (Clementini et al. 1999), 
and SN la abundance pattern by the W7 model. The Si 
abundances are clearly higher than the sum of the two. 

At the cluster center, the Fe/Si ratio and the Fe/O 
ratio of the gas of the Centaurus cluster are systemat- 
ically higher than those of M87, respectively. The main 
difference between M 87 and the Centaurus Cluster is 
the IMLR as described in the previous section. In figure 
M (Fe/S 

^)sNia ^ also plotted against IMLR. The derived 
(Fe/Si) SNIa value is much smaller than the value of 2.6 
in units of the solar ratio which is predicted by the W7 
model. Regions with larger IMLR have larger values of 
(Fe/Si) SNIa , when adopting the SN II pattern of either 
the Galactic stars (Clementini et al. 1999), or the nucle- 
osynthesis model (Iwamoto et al. 1999), A larger IMLR 
corresponds to a longer accumulation time scale. For M 
87, it needs only a few Gyr to accumulate the observed 
gas mass within lOkpc where the gas from the cD galaxy 
is dominant (Matsushita et al. 2002). In contrast, in the 
center of the Centaurus cluster, ejecta of SN la are accu- 
mulated over a much longer time scale than that of M 87. 
Therefore, the correlation between (Fe/Si) SNIa and IMLR 
support the suggestion by Finoguenov et al (2002) and 
Matsushita et al. (2003) from the M 87 observation that 
(Fe/Si) SNIa depends on the age of the system. Matsushita 
et al. (2003) discovered that the Si/O ratio of the Galaxy 
also increases systematically toward metal rich stars. This 
result indicates that there is a source of high Si/O ratio 
in our Galaxy at present and supports the scenario sug- 
gesting a dependency of the Fe/Si ratio of ejecta of SN la 
on the age of the system. 



stellar B-band luminosity of the cD galaxy (IMLR). The 
total B-band luminosity and the effective radius were cal- 
culated using Prugnicl and Heraudeau (1998), and the ef- 
fect of the Galactic extinction was corrected using Schlcgcl 
et al. (1998). The IMLR of the Centaurus cluster is much 
higher than that of the M 87, especially within the central 
region. Since the contribution of the SN la to Fe in the 
ICM is dominant, this indicates that accumulation time 
scale of SN la is higher in the Centaurus cluster than that 
in M 87, although we do not know the actual gas flow 
rate due to the cooling flow problem. The longer accu- 
mulation time scale in the Centaurus cluster means that 
M|n0sn/«* w as higher in the past. In other words, the Fe 
mass synthesized by a SN la was higher in the past or the 
ratio of SN la rate to stellar mass loss rate was higher in 
the past as suggested in Renzini et al. (1993). 



5.5. The hydrogen column density and the cooling 
flow model 

In the ASCA era, in order to explain observed fluxes of 
soft energy photons which are much smaller than the ex- 
pected value from the standard cooling flow model, ab- 
sorption within cluster cores was discussed (Allen et al. 
2001). In order to obscure a cooling flow in the Centaurus 
cluster from the standard model, we need at least ad- 
ditional column density around 2 x 10 21 cm~ 2 . However, 
with the XMM-Newton observations, the observed hydro- 
gen column density is too low to conceal the cool compo- 
nent from the cooling flow model (Matsushita et al. 2002, 
Bohringer et al. 2002). The observed hydrogen column 
density of the Centaurus cluster is also consistent with 
the Galactic value, 8.06 x 10 20 cm~ 2 over the whole field 
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Fig. 19. [Fe/Si] of the ICM of the Centaurus cluster (dia- 
monds) and M 87 (crosses; Matsushita et al. 2003) are 
plotted against [Fe/O]. The average value of Galactic 
metal poor stars (Clementini et al. 1999; open square) and 
the abundance ratio of SN II model using the nucleosyn- 
thesis model (asterisk) derived in Nomoto et al. (1997), 
assuming Salpeter's IMF (Iwamoto et al. 1999) are shown. 
The solid line and dashed line represent the relation of the 
abundance pattern synthesized by SN la with Fe/Si=l.l 
(the best fit relation of M 87) and Fe/Si=2.6 (W7 ratio; 
Nomoto et al. 1984), respectively. 

of view. This means that the soft energy photons from the 
Centaurus cluster are not obscured. 

6. Conclusion 

With the XMM-Newton observation, the abundance pat- 
tern of O, Mg, Si, and Fe of the ICM of the Centaurus 
cluster are derived. At the very center, the Mg and O 
abundances are close to the stellar metallicity. This result 
indicates that the gas at the very center is dominated by 
the supply of the stellar mass loss from the cD galaxy. 

The abundances of Si and Fe are close to each other in 
units of the solar abundance, while the O and Mg abun- 
dances arc smaller, like in the centers of other clusters 
and groups of galaxies. In order to explain the abun- 
dance pattern, the SN la ejecta from the cD galaxy should 
have a higher Si/Fe ratio than the standard W7 model by 
Nomoto et al. (1984). 

There are small differences in the Fe abundances and 
the Fe/Si ratio in the ICM between the Centaurus clus- 
ter and M 87 (Matsushita et al. 2003). These differences 
and the higher IMLR in the Centaurus cluster may re- 
flect a higher accumulation time scale of SN la. In order 
to confirm this suspicion, we will further study abundance 
patterns of the ICM of number of clusters and groups of 
galaxies with Suzaku which has a better energy resolution 
at the O lines. 
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